Single particle excitations of the Kondo lattice 
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We present a simple theory for the description of the single particle excitations of the Kondo 
lattice model. We derive an 'effective Hamiltonian' which describes the coherent propagation of 
single particle-like fluctuations on a strong coupling groundstate. Even for /-electrons which are 
replaced by Kondo-spins, the resulting spectral function obeys the Luttinger theorem including the 
/-electrons, and our calculation reproduces the complicated evolution of the spectral function with 
electron density seen in numerical studies. 
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The theoretical description of the Kondo lattice re- 
mains an unsolved problem of solid state physics, yet 
the solution of this problem is crucial for understanding 
the anomalous prperties of /-electron metals Q. The 
simplest model which incorporates the essential physics 
is the Kondo lattice model 
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Here we consider the 'minimal' model, where each unit 
cell contains two orbitals, one of them for the uncor- 
rected conduction electrons the other for the strongly 
correlated /-electrons. Then, cj (/J ) creates a con- 
duction electron (/-electron) in cell i, rii tCr =fJ ^fi^, and 
e(k)=J2i j exp(ik-(Ri — Rj)) t^j is the Fourier transform 
of the inter-cell hopping integral ti,j for c electrons. For 
simplicity we consider only the symmetric case, U = 2ef. 
For U^>V (|l|) can be reduced to 



H sc = ^2 e(fc)Cfe )CT Cfe !<T + S j)C • S 



(2) 
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where Si. c (<Sj,/) denotes the spin operator for conduc- 
tion electrons (/-electrons) in cell i and J = AV 2 jet. One 
problem which by many is believed to be at the heart of 
the solution is the way in which the more or less local- 
ized /-electrons, which are replaced by mere spin degrees 
of freedom in the strong coupling version (^) , partici- 
pate in the formation of the Fermi surface. De Haas-van 
Alphen experiments on heavy Fermion metals as well 
as computer simulations of Kondo lattice models ||^|] 
suggest that despite their 'frozen' charge degrees of free- 
dom the /-electrons participate in the Fermi surface vol- 
ume as if they were uncorrelated. The limiting cases 
V=0, J=0, which obviously do not allow for participa- 
tion of the / electrons in the Fermi surface, therefore rep- 
resent singular points, so that a perturbation expansion 
in the (small) parameters V or J may not be expected to 



give meaningful results. Rather, the interaction between 
/-spins and conduction electrons must be incorporated 
in a non-perturbative way, in a similar manner as the 
single-impurity Kondo effect || . It is the purpose of the 
present manuscript to present a minimum effort theory 
for the Kondo lattice which is based on this requirement 
and shows how the nominal participation of the localized 
electrons in the Fermi surface can be understood even 
in the complete absence of true hybridization. We de- 
scribe the system by an 'effective Hamiltonian' for the 
Fermion-like charge fluctuations on top of a strong cou- 
pling ground state, and show that this treatment leads 
to remarkable agreement with the scenario inferred from 
numerical simulations. 

As a starting point, we take the case of half-filling (i.e. 
two electrons/unit cell, corresponding to the 'Kondo in- 
sulator') and vanishing inter-cell hopping integrals tij. 
In this limit, the lattice ground state is simply the prod- 
uct of N single-cell ground states, each of them being 
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FIG. 1. Charge fluctuations and their propagation (left 
panel) and their representation in terms of 'model Fermions' 
(right panel). 

a two-electron singlet state (see the state (a) in Figure [l]) . 
In the following we consider this product state as a kind 
of 'vaccuum'. This state clearly is a total singlet, and has 
no magnetic order; it thus is appropriate to discuss the 
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paramagnetic phase of primary interest. Switching on 
the £y then produces 'charge fluctuations' on this vac- 
uum: an electron can jump from a cell i to another cell 
j, leaving the cell i in a single electron eigenstate with 
number v, and the cell j in a three electron eigenstate 
with number \i (see state (b) in Figure Q). We note that 
the Hamiltonian ([j]) allows for two such eigen states of a 
single cell, the strong coupling version (^|) only for one. 
In further steps, these charge fluctuations can propagate: 
an electron from the three-fold occupied cell j can hop 
to another neighbor j' leaving cell j with two electrons 
and j' in a three electron state (see Figure [I]c ) or, alter- 
natively, an electron can jump from another neighbor i' 
into cell i, leaving i in a two electron state, i' in a sin- 
gle electron state. When there is only one c orbital/unit 
cell, the single cell states of both 1 or 3 electrons are spin 
dublets, i.e. these states have the spin quantum numbers 
of ordinary electrons. As our basic approximation we now 
restrict the Hilbert space such that whenever there are 2 
electrons in one cell, they are in the singlet ground state. 
This restriction implies that the propagating charge fluc- 
tuations do not 'leave a trace' of excited cells, i.e. their 
motion under this constraint becomes completely coher- 
ent. We can then interpret the first step in Figure |l| as a 
pair creation process, where two 'book keeping Fermions' 
are created on nearest neighbors in the vacuum state, 
and the subsequent steps in Figure [l] as a propagation 
of these Fermions. More precisely, if cell number i is in 
the v th single electron eigenstate with z-spin a we in- 
terpret this as presence of a hole-like Fermion, created 
by aJ, CT , whereas the cell being in the /I th three electron 
eigenstate with z-spin a is modelled by the presence of 
an electron-like Fermion, created by bj^ a . Within our re- 
stricted Hilbert space the dynamics of the model then is 
decribed by the Hamiltonian H e f f=VH'P where 

h = J2(a <4,**i,* + a ) + E t ( v ^ h lA* 

+Vl, 3 a] a a^ + Vl^b) a b^) + H.c.]. (3) 

and V projects onto the subspace of states where no 
site is occupied by more than one Fermion. This kine- 
matic constraint reflects the fact that the state of a 
given cell must be unique, (for the sake of brevity we 
have suppressed the state indices /i, v in (§)). The pro- 
cedere is very much analogous to spin wave theory, with 
the sole difference that we are considering Fermionic 
charge fluctuations rather than Bosonic spin fluctua- 
tions. Due to the product nature of the basis states, 
the various matrix elements V, V and V" in (||) are 
easily expressed in terms of the electron removal and 
addition matrix elements of a single cell, mpES.v — 
(tt^ltvl*? 5 ) and m IPES ^ = (tf^ct |*$ 2) >. For exam- 
ple V(i,*/j,0»=*i,.j rnpES.u mi P ES,n ( nere ) denotes 
the n th eigenstate of a cell with n electrons). The 'on 
site energies' are defined as A v —E u — E . 



In analogy with linear spin wave theory we now relax 
the constraint enforced by V, whereupon the Hamilto- 
nian (^) is readily solved by Bogoliubov transformation, 

ll+^Y.v u k?\,*a„,-k,* + E P u k,i,a b l,k : v A Possibil- 
ity to treat the constraint in a more rigorous fashion 
would be to apply complete Gutzwiller projection to the 
wave functions obtained from (||); anticipating that the 
main effect is a renormalization of the matrix elements, 
we may expect that a calculation without the constraint 
at least qualitatively reproduces the correct one. Our 
main justification for this approximation, however, is the 
good agreement with exact cluster results as discussed 
below. As mentioned above, the Hamiltonian may 
be thought of as describing the coherent propagation of 
single particle-like fluctuations 'on top of the strong cou- 
pling 'vacuum' state. One may expect that there will also 
be e.g. spin-like fluctuations, which would correspond to 
a cell being occupied by two electrons in a triplet state. 
The creation of these spin excitations by the propagating 
charge fluctuations would be described by terms of the 
form a\ T Sj •o ! t ,t'%',t') where Sj is a bosonic spin-triplet 
operator. Such terms could be treated e.g. within the 
'rainbow diagram' approximation Q, but in the present 
manuscript we restrict ourselves to the coherent motion. 
For the strong coupling limit (||) the calculation becomes 
particularly simple: one finds mpES= m iPES= l/v 7 ^ 
whence V r i , i =-F/ j =^.=t i , i /2, A=A=3J/4, so that the 
dispersion relation reads 

E ± (k) = (1/2) [ e(k) ± VA 2 + e(fc) 2 ]. (4) 

Formally, this is equivalent to the hybridization of a dis- 
persionless 'effective' /-level in the band center with a 
free electron band with dispersion e(fc), the strength of 
the 'nominal' mixing clement being 3J/4. It should be 
noted, however, that the resulting energy gap of A does 
not arise from the formation of a bonding and antibond- 
ing combination of c-like and /-like Bloch states, as in the 
hybridization model; rather, this gap originates from the 
energy cost to break two intra-cell singlets in the course 
of a charge fluctuation. This gap therefore is of a similar 
nature as the energy gap in a superconductor. 
Having computed the eigenvalue spectrum, we proceed to 
the computation of the single particle spectral function. 
This requires to resolve the 'ordinary' electron operators 
in terms of the model Fermions. Taking into account our 
basic assumption, namely that a single cell with 2 elec- 
trons can only be in its ground state, we can expand e.g. 
the electron annihilation operator as 

<V = X/ m PES,vA,u,° + X, "'//V<,A;..- ■ ( 5 ) 

Then, Figure ^a shows the single particle spectral func- 
tion for a ID chain of the full Kondo lattice ([!]) at half- 
filling (see e.g. Ref. H for a precise definition of this 
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quantity) . Particle-hole symmetry implies that the Fermi 
energy is zero in this case. It is obvious at first sight, 
that our calculation is in excellent agreement with the 
results of exact diagonalization ||. To begin with, unlike 
any band theory approach, our calculation gives the cor- 
rect number of 'bands': the nearly dispersionless upper 
and lower Hubbard bands with almost pure / character, 
and the two 'hybridization bands' which resemble the re- 
sult for the strong coupling limit, (Q). The hybridization 
bands change their character from high intensity c-like 
to weak intensity /-like near 7r/2, i.e. the Fermi momen- 
tum of the unhybridized d electrons. Note that while the 
dispersion of these bands could be modelled by band the- 
ory, the sharp drop of spectral weight at ir/2 could not. 
In the Kondo limit, where the strength of the ' effective 
hybridization' is small, there are extended regions of flat 
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FIG. 2. Single particle spectral function for the ID Kondo 
lattice with e(k)=— 2 cos(fc), (7=10 and V= 0.5. The electron 
density is 2.0/unit cell in (a) (i.e. half-filling) and 1.6/unit cell 
in (b). Full lines (dashed lines) correspond to /-like (c-like) 
spectral weight. Peaks to the right (left) of the vertical line 
correspond to electron creation (annihilation). 

(i.e. 'heavy') bands. Quite obviously the basic idea of our 
approach, namely to broaden the ionization and affinity 
states of a single cell into bands works well as far as the 
dispersion of energy and spectral weight is concerned. 
The present approximation does not incorporate the in- 
frared divergences which occur in the impurity problem, 
and therefore does not reproduce the parameter depen- 
dences of the low energy scales || . For example the low- 
ering of the total energy/site as compared to the nonin- 
teracting case is ~ A 2 ln(W/ A) when calculated from (Q), 
and the /-like spectral weight Zf of the 'heavy bands' has 
the same functional dependence on the parameter values 
as in the single cell problem, i.e. Z~(V/e/) 2 . 
We proceed to the doped case and first consider the ex- 



pression for the electron number operator. We give ex- 
plicit expressions only for the strong coupling limit (||), 
but all considerations for the full Kondo lattice are com- 
pletely analogous. In the 'vacuum state', the number 
of electrons is 2N and the presence of an a-Fermion (b- 
Fermion) decreases (increases) the electron number by 1, 
so that the electron number operator should be simply 

N e = J^i a k,aa{ a + bl a b k ^ ) = ^ ll^alk,n,a- (6) 

This operator counts both localized and conduction elec- 
trons, so that reducing the electron density p e below 2 
will give a Fermi edge in the lower hybridization band 
which satisfies a 'nominal' Luttinger theorem, i.e. in- 
cluding the / -electrons. The physical origin, however, 
is the fact that we have a density 2 — p e of holes in the 
insulating 'singlet background'. This results in a 'hole 
pocket' of volume (1 — p e /2), which however is nominally 
equivalent to a Luttinger Fermi surface including the /- 
electrons. There is, however, an extra complication: the 
electron number should also equal the integrated PES 
weight, which is given by the expectation value of 

h.a 

The first term on the r.h.s. is the contribution of the 
/-electrons' lower Hubbard band. Using (^|), which now 
reads Cfe i<J =(l/v / 2)(a^. a + bk,a), this becomes 

N + 2 ^2( ak '° a k,a + h k,J>k,° + a k,* b l,a + b k,*ak,a)- 

fc,(T 

Counting the electrons in real space on one hand and inte- 
grating the spectral weight on the other hand thus results 
in different expressions for the electron number. This can 
hardly be a surprise, since in a strongly correlated elec- 
tron system spectral weight and band structure are de- 
coupled to a large extent. As an example let us consider 
the extreme Kondo limit, and assume that we are gradu- 
ally reducing the electron density from the 'Kondo insu- 
lator value' of 2. The chemical potential corresponding to 
the real-space electron count (^|) then will cut more and 
more into the 'heavy' band, to produce the nominal Lut- 
tinger Fermi surface volume. Since the spectral weight 
per fc-point is much smaller than 1 for the 'heavy' band, 
however, the amount of spectral weight which crosses 
from PES to IPES as the the Fermi momentum changes 
can never be consistent with the electron number. A 
simple rigid band picture therefore must fail to maintain 
the consistency of electron number and integrated spec- 
tral weight. To cope with this problem, we choose the 
simplest possible solution and enforce the consistency of 
'ordinary' electron count and spectral weight integration 
by adding both expressions, @ and (^), for the electron 
number to the Hamiltonian, each of them multiplied by a 
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separate Lagrange multiplier: H — > H + fiN e + XN' e . The 
notion of 'two chemical potentials' may seem awkward at 
first sight, but as we will show now, this approach results 
in a remarkable consistency with the numerical results. 
For the strong coupling limit (|2) this substitution gives 
the same dispersion relation as ( |) , however with the re- 
placements E(k)^E(k)—p, and e(fc)^e(fc)— A. The ef- 
fect of these modifications can be seen in Figure ^b, which 
shows the spectral function for the Kondo lattice away 
from half-filling. To begin with, the /-like upper and 
lower hubbard band remain unaffected by the change in 
electron concentration. The Lagrange multiplier \x acts 
like a standard chemical potential, which cuts into the 
'heavy' band and, as discussed above, produces a Fermi 
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FIG. 3. Momentum distribution for the strong coupling 
limit for different values of J, p c — 0.5. 

surface consistent with the nominal Luttinger theorem. 
On the other hand, the Lagrange multiplier for the spec- 
tral weight, A, modifies the dispersion relation far from 
the true Fermi momentum: the occupied part of the o 
like band structure shortens, the unoccupied part grows. 
In fact, the change of the c-like spectral weight alone is 
very reminiscent of what is expected for free (i.e. unhy- 
bridized) conduction electrons. This can be understood 
by carrying on the formal analogy of (|J) with a hybridiza- 
tion gap picture: A obviously plays the role of the 'on-site 
energy' of the effective / level. Next, in the Kondo limit 
J « f the momentum distribution for the conduction 
electrons, n c {k) has the limiting behaviour n c (k) w I 
for (e(fe) - A) < J, n c (k) w for (e(fe) - A) > J, and 
n c (k) = 1/2 for e(fe) = A. The constant energy surface 
e(fe) = A therefore represents a kind of 'pseudo Fermi 
surface', where n c {k) drops from nearly 1 to nearlyO over 
a distance of order A/vf (with vf the Fermi velocity of 
the conduction electrons). Since we have the sum rule 
(2/iV) J^k n c(k)—p c , the density of conduction electrons, 
it follows immediately that A ss E F , the Fermi energy for 
the unhybridizcd conduction electrons without counting 
the /-spins. In other words, the band structure of the 
Kondo lattice away from half-filling is equivalent to an 
effective /-level, which is pinned near the 'frozen core' 
Fermi energy for conduction electrons of density p c , and 
mixes with a matrix element of strength A. This non- 
rigid band like change of the spectral function upon dop- 



ing is precisely what was seen the cluster diagonalization 
results of Tsutsui et al. M. Next, Figure || shows the 
change of n c (k) with J for a ID chain of the strong cou- 
pling limit (||). For very large J (which is an unphysical 
limiting case) n c (k) is nearly constant and drops to at 
hp = 37r/4, the Fermi momentum for hybridized conduc- 
tion and / electrons. As J is reduced, the true Fermi 
surface discontinuity shrinks more and more, and the 
'pseudo Fermi surface' at k F =n/A, the 'frozen core' Fermi 
momentum, starts to develop. Again, this behaviour of 
n c (k) is in complete agreement with numerical results 

SB- 

In summary, we have outlined a simple approximate the- 
ory for the single particle excitations of the Kondo lat- 
tice model. It may be viewed as the construction of an 
effective Hamiltonian describing Fermionic fluctuations 
on a strong coupling vacuum state. While the calcula- 
tion requires a number of rather strong approximations, 
and consequently may fail to correctly reproduce the ex- 
treme low energy energy scales of the problem, we be- 
lieve that it also offers a number of advantages: despite 
its extraordinary simplicity the obtained results compare 
favourably with the secnario obtained in numerical simu- 
lations, and include some of the known constraints, such 
as the Luttinger Fermi surface and the paramagnetic sin- 
glet nature of the ground state; one may therefore hope 
hope that our approach captures the essential physics of 
the problem. The calculation moreover is independent 
of dimensionality or lattice geometry, and hence may be 
complementary to more sophisticated methods available 
in ID pi. Moreover, there are a number of obvious 
refinements, such as a more rigorous treatment of the 
kinematic constraint or the inclusion of spin-like fluctu- 
ations. 
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